Stationary Properties of a Randomly Driven Ising Ferromagnet 
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ON ' We consider the behavior of an Ising ferromagnet obeying the Glauber dynamics under the 

ON 1 influence of a fast switching, random external field. Analytic results for the stationary state are 

presented in mean-field approximation, exhibiting a novel type of first order phase transition related 
to dynamic freezing. Monte Carlo simulations performed on a quadratic lattice indicate that many 
^ , features of the mean field theory may survive the presence of fluctuations. 
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Many information processing systems, natural or artificial, have a macroscopic number of connected elements 
subject to external stimuli changing faster than the characteristic thermal relaxation time. However, attempts at 
O ' handling nonequilibrium stationary states of such systems have been made only recently Q[. As illustrated by the 
randomly driven Ising model presented in this Letter, stationary states far from equilibrium might be many times 
G more effective in dynamically storing information than thermal equilibrium states. 

The Ising ferromagnet in a time-dependent sinusoidally oscillating field has received recently a lot of attention, from 
both a theoretical and experimental point of view. On the theoretical side, Rao, Krishnamurthy, and Pandit (2| have 



presented a large A^-expansion of the cubic O(N) model in three dimensions and calculated the critical exponents 
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related to the area of the hysteresis loop. The underlying dynamic phase transition has been then studied within 
both mean-field || and Monte Carlo simulations jlf-Q] . The theory presented in this paper is a generalization of these 
ideas for the case when the external field is generated by chaotic dynamics and/or is a random variable. 

Let a spin system fl = (si, S2, . . . , Sj, . . . , sjv), Sj = ±1, be in local thermal equilibrium and denote by P(jl, t) the 
probability of finding the system in state fx at time t. The Master Equation used first by Glauber || for defining a 
stochastic dynamics for the Ising model reads then 

= J2 w(jt\jjk)P(pk;t) - P(m; t) J2 w{Ui\n) (1) 

> 

' where jli — (s%, S2, ■ ■ ■ , — s-j, ■ • ■ , %)• w(v\fl) denotes the transition rate from configuration jl into state v. Note that 
P(fi\ t) can be represented as a 2 N dimensional vector. It is useful to expand P(jl, t) in the orthonormal basis formed 

by all possible spin products, P(/x; t) = J2 a ir a Yl iea Si, where n a — ([\ iea Sj) t ||. In this 'spin-correlation basis' 
f***- ' the time-dependent distribution is given by n(t). The spectrum of the Liouville operator C 

On' 

- = ~C B(t) m (2) 

is invariant under such orthogonal transformations. The slowest relaxation time corresponds to the decay of the 
order parameter (the total magnetization) and is denoted by r sys . The typical phonon-spin interaction time is 
Tspin-fUp <S T sys . In what follows we assume that the external field B is a random variable sampled identically and 
independently at discrete times t n — tltb, (n = 0, 1, . . .) from the symmetric distribution: 



p(B(t n )) = \5{B{t n ) - B ) + \5{B{t n ) + B ) (3) 

As long as tb S> r sys the spin system relaxes to global thermal equilibrium. The situation is very different if 
T sys ^> tb Tsvin—flip- Now the stationary state is determined by the distribution of the external field. By 
integrating Eq. ([j|) over the (t n ,t n + tb) intervals, one obtains a 'coarse grained' discrete Master Equation 

K(tn+T B )-Tr(tn) ~ XU \ ( A\ 

= -LB(t n )K{t n ) (4) 

TB 

which still describes correctly the long-time behavior of Eq. (Q). Eqs. (^-^) map our problem into an iterated function 
system (IFS) @. 

The invariant probability density V s [^] induced by the dynamics (0) satisfies the Chapman-Kolmogorov equation: 
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V s {tt) = J dTf'V s {n') J dB p(B) 5(tt - c- £bTb n') = 5(n - &- LbTb ^) =KV s {^) (5) 

where [A] denotes the dynamic average J d,Tv'V s (ir') J dB p(B)A and JC the Frobenius-Perron operator. 

As an example we consider the mean field Ising model. In this case all multispin correlations factorize in the 
thermodynamic limit and the stationary distribution depends solely on the total magnetization. The energy is defined 
as usual, 

E = -j^^SiSj - n B B^Si (6) 

where /is is the Bohr-magneton. For the transition rate w(p,i\jl) we use the Glauber-form 

w(fii\P) = ^[1 ~ s. i tanh(— ^ Sj + H)] (7) 

where (3 = K — j3J, H — fipsB and a sets the time constant. Applying Eq. ([|) one obtains after performing 

the thermodynamic limit N — > oo: 

m(t + 1) = t&nh(Km(t) + H(t)), (8) 
where time is measured in units of tb- The field distribution Eq. (||) leads to the one-dimensional map 

!ta,nh(Km(t) + Hq) with prob. | 
(9) 
tanh(_RTm(t) — Hq) with prob. \ 

Since in the stationary state, Eq. (||), [m fc (t + 1)] = [m k (t)], using Eq. (JsJ) and simple algebraic manipulations we 
obtain that the fc-th moment of the stationary magnetization is given by 



\m k ] = 



v + h 



vh 



fc = l,2,... (10) 



where v — tanh(ii'm) and h = tanh(_ff). For example, when expanding up to fourth order in Hq = tanh(i?o) we get 
for the second moment 

|m2] = i-^(i-lg + 3 ^) <"» 

Higher moments can be calculated recursively in full analogy to methods introduced in Q . 

The two graphs below, Figs. ^ and ^, show how the map Eq. (^|) changes between high and low temperature. 
For further use let us denote by mi, 7712, and 7713 the possible fixed points of the equation m = tanh(Xm + Hq) in 
descending order. The stationary magnetization distribution undergoes a tangential bifurcation at the critical field 

H c = \\n—^-r+Km^ (12) 
2 1 + rrv 



where m2 — — = ± \J K K 1 for K > 1. The corresponding phase diagram is shown in the upper part of 
Fig. |3|. Below the phase transition the stationary magnetization distribution bifurcates into two symmetric, stable 
'spontaneous magnetization distributions' and a central repellor. The phase transition is first order, the average 
stationary magnetization jumps at the phase border. 

A different kind of transitions are related to the analytic structure of the invariant density. Following the notation 
introduced in p7| , one can identify a singular-continuous density with fractal support (SC-F) in both the paramagnetic 
and the ferromagnetic phase. When a gap opens between the upper and the lower branch of the map the invariant 
distribution has a fractal support with the capacity dimension do < 1. The border of the (SC-F) region is given by 
Krai = Hq in the para- and K(m\ + 7713) = 2Hq in the ferromagnetic phase. In the region between do = 1 and 
doo = 1_ the distribution is singular-continuous with Euclidean support (SC-E) j^j. Using the ideas developed in 
pfj| , we obtain d^ = 1 if K(l — ml) — \. The density distribution is absolutely continuous (AC) if all generalized 
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dimensions [ p"5| equal one, d q = 1, (g = 0, . . . , oo). These results are graphically summarized in the lower part of Fig. 
| 

The generalized free energy of such a driven system can be defined as — (3J- — A, where A is the largest Lyapunov 
exponent of the map Eq. (|^) . As expected, in the thermodynamic limit we obtain a dynamic average over the thermal 
mean field free energy at magnetization m. After performing the average over p(B) the generalized free energy per 
spin is given by 

~ W = J dmVs ^\ ln2(cosh(2A'm) + cosh(2H )) (13) 

Strictly speaking, Eq. ( |i~3| ) is the average free energy. When considering a finite system or a long but finite dynamic 
trajectory, the free energy is normally distributed. As shown in [jl2| for the one dimensional random field Ising 
model, in the SC-F region the multifractal spectrum can be directly related to the second cumulant of the free energy 
distribution. The arguments presented in [jl2j apply also to our case, a broad multifractal distribution leads to large 
free energy fluctuations. 

Consider now Fig. || at negative magnetization values. Close to but above the critical field (|lj) the upper branch 
of the map comes close but does not touch yet the diagonal. Let n be the average number of iterations needed to 
pass through this region starting from the lower corner by following the upper branch of the map. According to the 
general theory of type-I intermittency in one-dimensional maps, n ~ (Ho — Hc)^ 1 / 2 . Since in each iteration step the 
probability of jumping back along the lower branch is 1/2, the total time spent in the lower corner is proportional to 
2™, or 

T ~ 2" ~ 2 clH °- H - r ^ (14) 

with c ~ O(l) a constant. This behavior suggests a dynamic freezing transition (for more details see flj|| ). 

In order to test the predictions of the mean field theory in a more realistic setting, we performed Monte Carlo 
simulations for an Ising model with nearest neighbor interactions on a square lattice. The driving field is sampled 
from the distribution Eq. (||) after each Monte Carlo step (MCS). The left side of Fig. || shows the measured 
magnetization distribution in the paramagnetic phase, which is similar to Fig. H Below the critical temperature one 
obtains distributions similar to the ones displayed on the right side of Fig. |, to be compared with Fig. |[ Note 
that the critical field H C (K) is not a universal quantity and the square lattice values are different from the mean 
field ones. More details will be published elsewhere. Thermal fluctuations and finite size effects wash out the fine 
structure of the multifractal magnetization distribution predicted by the mean-field theory. However, the sharp peaks 
and the presence of gaps indicate that at least the main features of the magnetization distribution are preserved in 
two-dimensions . 

Besides the theoretical interest in describing such systems, we believe that our predictions can be tested with recently 
developed experimental techniques. Dynamic magnetization measurements have been recently performed in ultrathin 
Au(lll)/Cu(0001)/Au(lll) sandwiches or epitaxial Co/Au(lll) films [p^|-p0[. Similarly, hysteresis measurements 
on the ultrathin film Co/Au(001) |2lj] indicate that below T c these systems undergo a dynamic phase transition 
belonging to the Ising-universality class. More relevant to our theory, the time evolution of magnetization clusters 
can be optically recorded. The typical relaxation times range from minutes to a few seconds with increasing field 
amplitudes p0[ . This relatively slow relaxation rate allows for a simple experimental realization of the randomly 
driven external field. 

Ultrathin films are potential candidates for magneto-optical storage devices. At well chosen control parameters 
the stationary magnetization distribution of the RDIM displays several well separated peaks. Hence, by coding 
appropriately the time-sequence of field switches, one can - in principle - store locally more than two binary states. 

We are grateful to U. Ramacher and the ZFE, Siemens AG for the SYNAPSE1/N110 neurocomputer, on which the 
Monte Carlo simulations were performed. This work was partly supported by the DFG through SFB 517. 
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FIG. 1. Mean field map and the stationary distribution in the paramagnetic phase (K — 0.4 and Ho/K — 0.21). 
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FIG. 2. Same as in Fig. 1 but close to the critical field value (K = 2.0 and Ho/K = 0.266). Two disjoint distributions are 
created around the stable fixed points, a repellor in the middle. 
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FIG. 3. Mean field phase diagram. The upper part (Ho > 0) shows the border between the para- and ferromagnetic phase. 
In the lower part (Ho < 0) the regions denoted by SC-F, SC-E correspond to a singular- continuous invariant density with 
fractal and Euclidean support, respectively, while in the AC-region the density is absolutely-continuous. Note that the diagram 
is actually symmetric in Ho ■ 
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FIG. 4. Left: Magnetization distribution for the square lattice RDIM averaged over eight different initial conditions. 
K/K c = 0.4, K c = 0.4407, H /K = 0.5, lattice size 415 x 415. The simulation was run for 2 • 10 s MCS. Compare to 
Fig. 1. Right: Same parameters but K = 2K C , Hq/K = 1. Compare to Fig. 2. 
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